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Abstract. It is shown that linear time-dependent invariants for arbitrary multi- 
dimensional quadratic systems can be obtained from the Lagrangian and Hamil- 
tonian formulation procedures by considering a variation of coordinates and mo- 
menta that follows the classical trajectory and defines a noetherian symmetry 
transformation. 
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There are several physical systems which in different approximations can be rep- 
resented by Hamiltonians which are quadratic forms of coordinates and momenta 
variables, i.e., the classical small oscillations problem. Thus the applicability of 
multidimensional quadratic Hamiltonians ranges from solid-state physics to particle 
physics, as for example in the relativistic oscillator models of elementary particles 
[1,2] and the particle creation in non-stationary metrics [3]. The exact solutions for 
stationary multidimensional quadratic systems have been found and discussed in the 
framework of canonical transformations [4]. 

Recently for the classical small oscillations problem by means of non-noetherian 
transformations, noetherian quadratic time-independent integrals of motion have 
been found [5]. For the one- mode parametric classical oscillator Ermakov found 
an invariant, which is quadratic in the momentum and position variables [6]. In the 
last years on the base of a geometric approach related to the Noether's theorem the 
last invariant has been rederived [7]. 

In the case of non-stationary quantum systems there are different procedures for 
solving time-dependent Hamiltonians and calculating evolution operators, but the 
simplest, at least for quadratic systems, is the method of time-dependent integrals of 
motion [8]. 

The theory of integrals of motion which do not depend on time in Schrodinger 
representation is well known [9]. On the other hand, for the one-dimensional quantum 
oscillator a time-dependent integral of motion quadratic in position and momentum 
operators has been found [10], which in form coincides with the Ermakov invariant. 
The linear time-dependent invariants have been found and applied to get the solutions 
for the quantum parametric oscillator and a charge moving under a time-dependent 
magnetic field [11]. In the most general multidimentional case, the time-dependent 
invariants of the non-stationary quadratic quantum systems also have been found 
and used to get the solutions of the time-dependent Schrodinger equation in the 
form of gaussian coherent states and Hermite polynomials of several variables [8, 12]. 
These invariants are linear combinations of coordinates and momenta operators, with 
time-dependent coefficients. 

We think it is necessary to understand the existence of these time-dependent 
constants of motion within the framework of a variational formulation. Up to now 
for one- mode parametric oscillator, the linear invariants of Ref. [11], were obtained in 
[13], [14] and [15] through the Noether's theorem procedure. Also with this formalism 
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in Refs. [13] and [14] the quadratic time-dependent invariants of this system were 
derived. Nevertheless until now it was unknown what symmetry corresponds to linear 
time-dependent integrals of motion in the frame of Noether's theorem prescription 
for multidimentional parametric oscillators. 

The aim of this paper is to get the linear multidimensional time-dependent in- 
variants by means of the Lagrangian and Hamiltonian formalisms of the Noether's 
theorem procedure. We give explicitly such a variation or noetherian symmetry 
transformation, which has a clear physical interpretation. Preliminary results, of the 
Lagrangian Noether's theorem procedure, for an arbitrary mutidimensional forced 
parametric oscillators with a detailed discussion of a two dimensional quantum sys- 
tem was presented in Ref. [15]. 

Next we present the active viewpoint of the Lagrangian and Hamiltonian for- 
malisms of the Noether's theorem procedure for a given Hamiltonian system. 

Lagrangian Formulation 

Let us consider an arbitrary time-dependent Hamiltonian, which is given in terms 
of a general quadratic form, i.e., 

H = ^Q a B ab (t)Q b + C a (t)Q a , (1) 

where hereafter latin indices run from 1 to 2n and sum over repeated indices is 
understood. Besides, we defined the coordinates and momenta like a 2n-dimensional 
vector 

/Pi\ 



Q 



Pn 
Ql 



(2) 



The 2n-dimensional matrix B and vector C, which define the quadratic form of the 
Hamiltonian, can be rewritten in terms of the n x n submatrices, A, B, C, and D, 
and the n— column vectors F, and G, respectively 

The hermiticity of the Hamiltonian implies that the matrix B must be symmetric, 
and this gives rise to the following symmetry conditions over the four constituents 
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n x n matrices: 

A 1 = A, B l = C, D l = D . (4) 
Expanding the Hamiltonian (1) we obtain 

H = - (A a pp a pp + 2B a pp a qp + D a pq a qp) + F a p a + Gpqp , (5) 

where the symmetry conditions (4) were used. Now and hereafter, greek indices run 
from 1 to n. From Hamiltonian (5), it is straightforward to obtain the Hamilton 
equations of motion 



Pa = -CapP/3 - Dapqp - G a , 
4a = A^ppp + Bapqp + F a . 



(6a) 
(66) 



This coupled set of first order differential equations can be decoupled in the following 
sets of second order equations 



Pa = ~ 
+ 



C - ( D — DB | D~ x 



Pf3 ~ 


C + DA- ( 


. a/3 





P(3 



a/3 



q a = 



D + DB^jD- 1 



B + (A- AC) A' 1 



G/3 + D a/3 F{3 + G a , 



(7a) 
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B-AD-{ A- AC^jA^B 



a/3 



A- AC 



A- 1 ] F/3 - A a pGp + F a 

a/3 



(7b) 



Next, we are going to find the Lagrangian of the system and apply Noether's 
procedure to get the time-dependent invariants. The Lagrangian is obtained by 
making a Legendre transformation to the Hamiltonian (5) and using the equation of 
motion (6b), resulting 



L =\ {A-\ p q a qp — (D — C A~ l B) a pq a qp + A^F^} 
- (A- 1 B) af3 q a q f3 - A~\F a qp - G a q a + (A~ 1 B) a pF a qp 



(8) 



To get the constants of the motion for this system, let us propose infinitesimal vari- 
ations of coordinates given by 

Sq a = h a (t) , (9) 



where h(t) is an arbitrary n— dimensional vector depending only on time. The cor- 
responding variation induced in the Lagrangian (8) is 

SL =A~ph a (qp - Bp y q 7 - Fp) - (CA'^aphaqp 

+ {CA~ X B - D) af3 h a qp - {CA- 1 )^^ - h a G a . 

This expression can be rewritten as 

6L=± 
dt 



KlKq, - (CA-%^ - f * (A-lKF, + (CA-% M + K G a ) } 
(A-*h a y + (A-^pha - {(CA-XpK)- - (CA^B - D) a ph a ] 



qp ■ 
(ii) 

Asking that the variations of coordinates (9) yield a Noether's symmetry transfor- 
mation, one has that (11) must be equal to a total time derivative of a function 
^L(q a ,t). This is achieved if the variation h(t) satisfies the differential equation 

(h a A-iy + h a (A- 1 B) af3 
- (KiCA- 1 )^)- - h a {CA- x B - D) a p = , 

and the function 0,L(q a ,t) is defined by the curly brackets of expression (11). De- 
veloping (12), one sees that it has the same form that the homogeneous part of the 
differential equation of motion for q a , given in (7b). For these symmetry transforma- 
tions according to the Noether's theorem the associated conserved charges are given 
by Jl = (dL/dq a )5q a - fi L , that is 

Jl = Pah a - A'phaqp + (CA'^aphaqp 

+ J dtfA-lhaFp-iCA-^aphaFp + haGa} , (13) 

with p a = [A~^qp — (A~ 1 B) a pqp — A~pFpJ. The different conserved charges are 
obtained by the initial conditions estabished for the system (12). 

Hamiltonian Formulation 

There is also the possibility of using the Noether's theorem in the Hamiltonian 
formulation. In this case we consider the action 

S = I' dt h: +ab Q b Q a - H(Q a , t)\ , (14) 

Jto 

5 



where we have introduced the (2n x 2n)-matrix in block form: 

E+= (o o) ' 

with I being the identity matrix in n dimensions. 

Taking an arbitrary variation of (14) with respect to coordinates and momenta 
defined by (2) we get the expression 

SS = ^ 1 dt joQ a (v +ab Q b - + 5Q a E_ ab Q 6 | , (15) 

where the transpose of the matrix E + was used, 



£_ = 
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If the transformation 5Q a is a noetherian symmetry , the curly bracket in (15) must 
be written as a total time derivative of a function 0#(Q a , t), i.e., 

Z+abQb - -qqA + SQa^-abQb = ^^#(<2a, t) . (16) 

Following now the same procedure as in the Lagrangian case the associated conserved 
quantity to this variation is given by 

Jh = 5Q a T l - ab Q b - fl H ■ (17) 

This can be checked also by taking the total time derivative of (17) and the Hamilton 
equations of motion together with the condition (16). 

For the multidimensional quadratic parametric Hamiltonian (1), the linear time- 
dependent constants of motion can be obtained by means of the variations 

SQ a = (Sp a , Sq a ) = (g a (t) , h a (t)) . (18) 

Substituting these variations into the left hand side of (16), we find a symmetry 
transformation if the following relations are satisfied: 

g a = -C a pgp - D aP hp , (19a) 
ha = A a pgp + B a php , (196) 
&h = g a q a - J dt{F a g ol +G ol h ol ) . (19c) 
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It is important to remark that the first two equations have the same form that the 
homogeneous part of the Hamilton equations of motion given in (6). Therefore it is 
straightforward to demonstrate that the variations h a satisfy the same differential 
equations (12) as the corresponding variations of the Lagrangian formulation. For 
this reason, the constant of motion (17) coincides with its Lagrangian counterpart 
(13), i.e., J L = Jh- 

The expression (13) gives rise to 2n invariants correponding to the 2n indepen- 
dent solutions for the system of equations (12). These 2n constants of the motion 
can be rewritten in matrix form 



The constants of motion defined in (20) satisfy at t = the conditions po(0) = p and 
<fo (0) = q. These requeriments imply different initial conditions for the system of 
differential equations (12), and guarantee the existence of 2n independent solutions, 
which are denoted by the superscript on the function h(t). Besides these vector 
solutions are written horizontally. The time-dependent column vector A is given by 



It is important to remark that the expressions in Eq. (20) are the linear invariants 
proposed in Ref. [8,12], where it was shown that A is a 2n dimensional symplectic 
matrix, as it must be because the constants of motion po(t) and qo(t) satisfy the 
canonical commutation relations. However now the elements of the matrix A are 
given in terms of the 2n independent solutions of the homogeneous classical equations 
of motion (12), and their derivatives. 

To have that the invariants coincide with the position and momentum operators 
at t = 0, the appropriate initial conditions for the solutions of Eq. (12) are 




(20) 



where the transformation matrix A takes the form 




(21) 




(22) 




1 < i, j < n 

n + 1 < i, j < 2n 



(23a) 
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and for their derivatives 



1 < i, j < n 

n + 1 < i, j < 2n 



(236) 



So we have proved that for multidimensional stationary or non-stationary classical 
or quantum systems the linear time-dependent invariants follow from the Lagrangian 
or Hamiltonian formulations of the Noether's theorem procedure by considering a 
variation of coordinates or momenta that follows the classical trajectory. 

An additional remark is that the quadratic invariants for the stationary clas- 
sical small oscillation system found from non-noetherian transformations [5] can be 
obtained from the linear invariants presented here. 

The suggested method can be extended to classical or quantum field theory by 
considering infinite degrees of freedom and as consequence finding an infinite num- 
ber of time-dependent invariants linear in fields and their corresponding conjugated 
momenta. 

One of the authors (V.I.M.) would like to thank Instituto de Ciencias Nucleares, 
UNAM for its hospitality. 
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